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Abstract 



The theory of the tight span, a cell complex that can be associated 
to every metric D, offers a unifying view on existing approaches for 
analyzing distance data, in particular for decomposing a metric D into 
a sum of simpler metrics as well as for representing it by certain specific 
edge- weighted graphs, often referred to as realizations of D. Many of 
these approaches involve the explicit or implicit computation of the 
so-called cutpoints of (the tight span of) D, such as the algorithm for 
computing the "building blocks" of optimal realizations of D recently 
presented by A. Hertz and S.Varone. The main result of this paper is 
an algorithm for computing the set of these cutpoints for a metric D 
on a finite set with n elements in 0(n 3 ) time. As a direct consequence, 
this improves the run time of the aforementioned 0(n 6 )-algorithm by 
Hertz and Varone by "three orders of magnitude" . 

Keywords: metric, cutpoint, realization, tight span, decomposition, block 

1 Introduction 

The decomposition of a given metric into simpler metrics (see e.g. [I5j) is a 
fundamental problem in classification featuring applications in, for example, 
clustering (e.g. [2]), "networking" (e.g. [3]), and phylogenetics (e.g. [13]). 
The theory of the tight span 

T(D) := {feR x : f(x) = sup(D(x,y) - f(y)) for all x G X}, 

y ex 

of a metric D denned on a set X [T71 [B] offers a unifying view on various 
existing approaches developed for this task. In this paper, we focus on de- 
compositions of metrics D defined on a finite set X that are induced by 
cutpoints of T(D), that is, maps / G T(D) such that T(D) — {/} is discon- 
nected. These decompositions are closely related to certain graph realizations 
of D, that is, connected edge-weighted graphs G = (V, E, £ : E — > M>o) with 
X C V for which D(x,y) = Dc(x,y) holds for all x,y G X (where Dq 
denotes the shortest-path metric induced by G on F). 

To describe these graph realizations, recall (see e.g. [20]) that a vertex v 
in a graph G = (V, E) is called a cut vertex (of G) if there exist vertices 
u,u' G V with {u,v}, {u',v} G E such that every path from u to u' in G 
passes through v. Moreover, a maximal subset Bey with the property that 
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Figure 1: (a) An example of a metric D on X = {a, . . . ,e}. (b) A block 
realization of D: The vertices in the shaded region form a block and edge s 
is a bridge, (c) A block realization of the restriction D' of D to the subset 
X' :=X\{c}. 



the induced graph G[B] := (B,E fl ( 2 )) has no cut vertex is called a 
of G. A graph realization G of D is called a realization of Z) if G is a 
graph, i.e., every block of G is a clique, and every vertex in V \ X has 
degree at least 3 and is a cut vertex of G. An example of a block realization 
is presented in Figure [T^b) . 

In a recent series of papers [HI EJ [TO] , it has been observed that a map / G 
T(D) is a cutpoint if and only if the graph Tf : = (Xf, Ef) defined, for every 
/ G R x , by X f := supp(/) and E f := {{x,y} G ( sup 2 p(/) ) : f(x) + f(y) > 
D(x,y)} is disconnected (where, as usual, supp(J) := {x G X : /(x) 7^ 0} 
denotes the support of /), that a map / in 

P(D) := {/ G M' Y : f{x) + /(y) > D(x, y) for all a; G X} 

for which the graph is disconnected must be contained in — and, hence, 
a cutpoint of — T(D), and that a cutpoint / G T(D) has a neighborhood 
that is homeomorphic to the open interval (— 1,+1) if and only if Tf is 
the disjoint union of two cliques. As such maps are of little interest for 
constructing block realizations, we will focus our attention in particular to 
the set of those cutpoints, denoted by cut*(D), for which either supp(J) 7^ X 
holds or Tf is not the disjoint union of two cliques. 

In this paper, we present an algorithm with run time 0(n 3 ) to compute 
cut*(D), where n = \X\, improving the run time of the algorithm presented 
in [7j. Once the set cut*{D) is available, it is straight-forward to compute 
a corresponding canonical block realization G = Gd = (Vr>, Ed, £d) of D in 
0(n 3 ) time. And, from that block realization it is then easy to derive, for 
every block B of Go, a corresponding metric D B on X that assigns, to any 
two elements x,x' G X, the value Db{x,x') defined as the total weight of 
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those edges on any shortest path from x to x' in G that have both end points 
in B. For example in Figure QJb) the distance DB{a,d) between a and d is 
5 where B is the block in the shaded region. 

This yields a decomposition of D into a sum of metrics of the form Db 
where B runs through the blocks of Gd that can be computed in 0(n 3 ) time. 
As a consequence, our algorithm improves the run time of the algorithm 
presented in [T4] that follows a 2-step approach: It computes first those 
metrics Db that correspond to blocks B with only 2 vertices, the so-called 
bridges, (see [13] for details) and then the remaining metrics Db for the 
blocks B that are not bridges. 

Our paper is structured as follows. In the next section, we introduce 
the concept of block splits and show how they can help in the computation 
of cut*(D). In Section [3], we establish the key properties of block splits and 
cutpoints in cut*(D) that we use in our new algorithm for computing cut*(D), 
and we present this algorithm in Section HI 

2 Block splits 

In this section, we present a key concept used in our algorithm for computing 
the set cut*(D), where D is the given metric on a finite set X: Recall that 
a split S of X is a bipartition of X into two non-empty subsets A and B, 
also denoted by A\B or B\A. A split A\B of X is called a block split of 
X (relative to D) if there exists a map / G P(D) with supp(/) = X such 
that is the disjoint union of two cliques whose vertex sets are A and B, 
respectively. Note that the condition used in the definition of a block split 
above is slightly stronger than the condition given in [TBI p. 10]. Also note 
that block splits, although not given a specific name, play an important role 
in the algorithm for computing bridges presented in [13]. The set of block 
splits of X induced by D is denoted by S^. In the following, we will also 
often simply write xy for D(x, y), x, y G X. 

Our first goal is to establish a property of block splits that allows to 
efficiently check whether a given split of X is a block split. To this end, we 
first recall the following well-known fact. 

Lemma 2.1 Given two sets A and B and a bi-variate map (p '■ A x B — ► R 
from the Cartesian product A x B into the real numbers (or any Abelian 
group), there exist maps 4>a '■ A — > R and 4>b '■ B — »■ M with cf)(a, b) = 4>a{q) + 
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4>B{b) for all a G A and b G B if and only if 0(a, b)+<p(a' , b') = (p(a, b')-\-(j){a! ', b) 
holds for all a, a' G A and 6,6' G I? if and only if 4>(a,b) + 0(ao,&o) — 
0(a, 60) + 0(^0) holds for some fixed elements ao G A and bo & B and all 
a G A and b E B. 

Proof: If there exist maps 0^ : A — > K. and 0£ : -B — > M with 4>(a,b) = 
</>a(&) + 4>B{b), one clearly has 0(a, 6) + 6') = </>A( a ) + + 0a(°O + 

4>B(b') = 4>(a, b') + 0(a', 6) for all a, a' & A and b,b' E B while, conversely, if 
4>(a, b) + 4>(ao, b ) = <p(a, b ) + 0(ao, 6) holds for some fixed elements a E A 
and bo £ B and all a G A and 6 G B, one has (f>(a,b) = 0a(o) + 4>B{b) for, 
e.g., the two maps 0^ : A — ► R : a 1— ► 0(a, &o) and 0# : -B — > K. : & 1— > 
0(oo,6) - 0(oo,6o)- ■ 

Next, we define, for any map / G P{D) and any subset Y of X, the 
virtual distance D(f\Y) from / to Y (relative to D) by 

D(f\Y) := i min{/(y) + f(y') - yy' :y,y'EY}. 

We will also write D(x\Y) rather that D(f\Y) in case / coincides with the 
so-called Kuratowski map k x [18] associated with an element x E X defined 
by k x (y) := xy for all y G X. Note that < D(f\Y) holds for all / and Y as 
above. Note also that, given a split S = A\B of X with ab + a'b' = ab' + a'b 
for all a, a' G A and b, b' G B, and any two elements a G A and b £ B, one 
has 



D(a|5) + D(6|A) - a& (1) 
= - min {ab' + ab" + a'b + a"b - b'b" - a' a" - 2ab] 

2 a',a"£A;b>,b"eB 

= - min {a'b' + a"b" - b'b" - a' a"} 

2 a',a"£A;b',b"eB 

= - min {max(a'&' + a"b", a'b" + a"b') - a' a" - b'b"} =: a s , 

2 a',a"€A;b',b"€B l K ' ' 



and that as has been dubbed the isolation index of S pQ. 

To illustrate the above definitions, note that, for the metric given in 
Figure (Ha), the split S = {a, b}\{c, d, e} is a block split with D(a\{c, d, e}) = 
3, £>(&|{c,d,e}) = 2, D(c|{a,6}) = 4, D(d\{a,b}) = 7, and D(e|{a,b}) = 3 
and, therefore, D(x\{c, d, e}) + £)(?/|{a, 6}) — D(x,y) = 1, for all x G {a, 6} 
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and y G {c, d, e}, the weight of the edge s separating {a, b} from {c, d, e] in 
the corresponding block realization depicted in Figure [CJb). 
More generally, we have 

Lemma 2.2 A split S = A\B of X is a block split of X if and only if the 
isolation index as of S is positive and, choosing arbitrary elements a® G A 
and bo G -B, ao&o + = &o&' + a 'bo holds for all a! E A and 6' G -B. 

Proof: Assume first that S = A\B is a block split. By the definition of a 
block split, there exists a map / G P(D) for which is the disjoint union 
of two cliques whose vertex sets are A and B and, therefore, we clearly have 
D(f\A), D(f\B) > 0. Moreover, for the restrictions 4>a '■= /U an d 4>b '■= f\s 
of / to A and B, respectively, we have 4>a(o) + 0b(&) — ao for all a G A and 
b £ B, and, therefore, afr + a'b' = ab' + a'b must hold for all a, a' £ A and 
b,b' E B in view of Lemma [2.11 applied to the bivariate map := D\axb- In 
consequence, by Equation ([1]), we have D(a|S) + D(b\A) — ab = as for all 
a <E A, b <E B and, so, choosing any a G A and b E B, we also have 

a 5 = L>(a|5) + £>(&|A) - a 6 

= \ , „ mm {/(&') + f(b") - b'b" + f(a') + f(a") - a'a"} 

Z a,a"£A;b,b"£B 

= D(f\B) + D(f\A) >0 , as required. 

Conversely, choosing arbitrary elements a G ^4 and 6 £ B, if a 6 +a / 6 / = 
aofr' + a'bo holds for all a' E A and b' E B and the isolation index of 5 is 
positive, then, in view of Lemma 12.11 we may choose any two non-negative 
real numbers ^a, 1b with ja + 1b = «s and define the map 

/ = fA->-Y A ,B-+~f B ■ X — ► R (2) 

by f(a) : = D(a|S)- 7A for all a G A and f(b) := D(b\A)- lB for all b G 5. In 
view of Equation (Tj[|), we clearly have f(d)+f(b) = D(a\B)+D(b\A) — a s = ab 
for all a G A and b £ B. Moreover, we have f(a) + /(a') > aa' for all 
a, a' G A as /(a) + f(a') > D(a\B) + D(a'|B) - 2a s = ab + a'b - 2D(b\A) = 
aa' + (ab + a'b — aa') — 2D(b\A) > aa' holds for all a, a' G A and every b £ B 
in view of the definition of D(b\A) (indeed, ab + a'b — aa' is one of the terms 
whose minimum, over all a, a' G A, coincides with 2D(b\A)j, and we have 
f(a) + f(a') > aa' for all a, a' & A if and only 7^ < as holds as ja = &s 
implies f(a) + f(a') = aa' for all a, a' £ A with a& + a'b — aa' = 2D(6|A). By 



6 



symmetry, we have also f{b)+f{b') > bb' for all b,b' G B and fib) +f(U) > bb' 
if and only if 7# < as holds. Thus, Ef is a subset of (^) U (^) , and it coin- 
cides with this set if and only if 7a, 7s < cts holds. So, A\B must indeed be 
a block split, as required. ■ 

It is also worth noting that, for every block split S = A\B, every / G 
P(D) with E f C [fj U (f) (or, equivalently, with f(a) + /(6) = ab for 
all a G A and b E B) actually is of the form / = fA^ lA ,B^ B for some 
7a, 7b > with 7a + 7b = as: Indeed, in view of Equation ([I]), we have 
D{a\B) - f(a) = ab + a s - D(b\A) - f{a) = f(b) + a s - D{b\A) for all 
a G A and b G B in this case, implying in particular that neither side 
changes once we replace a by any other element in A nor b by any other 
element in B. So, choosing any fixed a G A and &o G -B, we may put 
7a := D(a \B) — f(ao) and 7b := D(b \A) — /(&o) in which case we have 
7a+7b = D(a \B)+D(b \A)-f(a )-f(b ) = D(a Q \B)+D(b \A)-a b Q = a s , 
f(a) = D(a\B) - 7a and f(b) = D(b\A) - 7 B for all a G A and b G B. 
Moreover, we have 7a > in view of D(ao\B) = ^minjao^ + a^b' — bb 1 : 
b,b' G B} = /(a ) + imin{/(6) + /(&') - 66' : 6, 6' G 5} > /(a ), and, 
similarly, 7s > 0. 

In other words, given any block split S = A\B, the set 



T(D\S) := {feP(D):E f C 




forms an straight line segment in IR parametrized by the straight line seg- 
ment {(7a, 7b) G IR>o : 7a + 7b = as 1 } in K 2 , and the two end points 
f A ■= f'A^a s ,B^o (closer to A) and f B := fA->Q,B-, as (closer to B) must 
each be either cut points of T(D) that do not have a neighborhood that is 
homeomorphic to the open interval (—1, +1) or elements of the set K(D) : = 
{k x : x G X} consisting of all Kuratowski maps for D. Hence we have the 
following. 

Corollary 2.3 For every block split S = A\B the maps /a and fs must be 
contained in the set Cut*(D) := cut*(D) U K(D). 

3 Key properties of and Cut*(D) 

As we have seen in the previous section, it is sometimes helpful to consider 
the bigger set Cut*(D) rather than cut*(D). Since we can easily identify 
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those Kuratowski maps that are not in cut*(D), we will now focus mainly 
on Cut*(D). The following lemma establishes the key properties of and 
Cut*(D) that we will use in our algorithm to compute these sets recursively. 

Lemma 3.1 Let x be an arbitrary element of X. Define X' :— X \ {x} and 
let D' denote the restriction of D to X'. Then the following assertions hold. 

(i) If S = A\B is a block split of X, then either S = {x}\X' or the 
restriction A n X'\B n X' of S to X' is a block split of X'. 

(ii) If / G Cut*(D)\K(D) = cut*{D)\K(D) has the property that there is 
no block split S = A\B of X with / G {/a, /b}, then the restriction /' of 
/ to X' is an element of Cut*(D') and /(x) = max {xy — f'(y) : y G X'} 
holds. 

Proof: (i) Clearly if S = A|5 is a block split of X with A,5 ^ {a;}, then 
5" = A fl fl I' is a split of X', and there exists a map / G P(L>) 

such that is the disjoint union of two cliques with vertex sets A and B 
implying that the restriction /' of / to X' is in P(D') and that Tji is the 
disjoint union of two cliques with vertex sets AnX' and BC\X', respectively. 
This establishes (i). 

To see that (ii) holds, suppose / G cut*(D) \ K(D) and that there is no 
block split S = A\B of X with / G {/a, Jb}- Clearly, the restriction /' of / 
to X' is in P(D'). So, it remains to show that ry is disconnected, but not 
the disjoint union of two cliques, which implies in particular that /' G T(D'). 

Assume for a contradiction that Tf> is connected or the disjoint union of 
two cliques. We first note that this implies that T f has at least one connected 
component that is a clique. To see this, observe that if V f> is connected, then 
Tf has precisely two connected components, one of whom consists of the 
single vertex x, thus trivially forming a clique. Similarly, if is the disjoint 
union of two cliques, then one of these cliques is also a connected component 

ofiy 

Let A denote the vertex set of a connected component of Tf that forms 
a clique. Note that this implies that D(f\A) > 0. Put B := X \ A and 
S := A\B. We want to show that S is a block split with / G {/a,/b}j 
yielding the required contradiction. To this end, choose arbitrary elements 
o G A and b G B. Since / G cut*(D), B cannot be the vertex set of a 
clique in Tf, and so there must exist two distinct elements bi,b 2 G B with 
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the property that fibx) + f(b 2 ) = bib 2 , implying that D(f\B) = holds. 
Since 

a'b + ab' = f(a') + f(b) + f(a) + f(b') = a'b' + ab 

clearly holds for all a, a' G A and b, b' G B, we have, in view of Equation (1TT) 
and the definition of Tf, 

a s = D(a \B) + D(b \A) - a b 

= (/(oo) + D(f\B)) + (/(& ) + A)) - /(oo) - /(bo) 
= £>(/|il)>0, 

and, therefore, S 1 is indeed a block split. 

It remains to show that / G {/a, Jb}- More specifically, we will show that 
f = Jb holds. By the definition of fs and in view of the fact that D(f\B) = 
and D(f\A) = a s holds, we have indeed f B (a) = D(a\B) = f(a) + D(f\B) = 
f(a), for every a G A, and f B (b) = D{b\A)-a s = f{b) + D{f\A)-a s = f{b), 
for every b G B, as claimed. ■ 

We close this section with establishing bounds on the size of the sets E p 
and Cut*(D) that we will use in the analysis of the run time of our algorithm 
in Section HI 

Lemma 3.2 Let D be a metric on a finite set X with n elements. Then 
|E D | < 2n - 3 and \Cut*(D)\ < An - 5 holds. 

Proof: To establish the first claim, it suffices to note that any two splits 
A]\Bi, A 2 \B 2 G T,£) are compatible, that is, at least one of the four intersec- 
tions Ai fl A 2 , Ai fl S2, Si fl A 2 and B 1 fl -B 2 is empty, since it is well known 
that every set of pairwise compatible splits of X contains at most 2n — 3 
splits (see e.g. Proposition 2.1.3 and Theorem 3.1.4 in [19]). So, assume 
for a contradiction that there exist two splits Ai|Si and A 2 \B 2 in that 
are not compatible. Then we can choose arbitrary elements a G A\ fl A 2 , 
b G Bi n A 2 , c G A x n B 2 and d G S x n B 2 . By the definition of a block split, 
there exist maps fa G T(D), i G {1,2}, for which the graph Y ^ is the disjoint 
union of two cliques with vertex sets Ai and Bi. But then, by the definition 
of T f ± and T f 2 , 

fx(a) + fx(b) + fx(c) + fx(d) = ab + cd< f 2 (a) + f 2 (b) + f 2 (c) + f 2 (d) 

= ac + bd< fx(a) + fx(b) + fx(c) + fx(d) 
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holds, a contradiction. 

Next we show \Cut*(D)\ < 4n — 5. Since, clearly, \K(D)\ < n, it suffices 
to show that \Cut*(D) \ K{D)\ < 3n — 5. In [8], it is shown that there 
exists a block realization G = Go of D such that the cut vertices in G are in 
one-to-one correspondence with the elements in Cut*(D) \ K(D). Moreover, 
the number of cut vertices in any graph is well known to be less than the 
number of blocks of this graph (see e.g. [12] )■ Hence, it suffices to show that 
the number of blocks in G is at most 3n — 5. Yet, it has been shown in [9] 
that there is a canonical bijection from the set of blocks of G to a set II of 
(strongly) compatible partitions of X, that is, of partitions such that there 
exist, for any two distinct partitions %i and 7r 2 , two necessarily unique subsets 
A\ G tti and A 2 G 7r 2 of X with A\ U A 2 = X (generalizing the concept of 
compatibility for splits to arbitrary partitions of X). Therefore, it suffices 
to show that, for all n > 2, every set of pairwise compatible partitions of X 
contains at most 3n — 5 partitions which we will establish by induction on 
the size of X. Clearly, if n — 2 then there is only one partition of X. 

Now assume n — \X\ > 2. If every partition in II is a split of X, then 
|n| < 2n — 3 < 3n — 5 must hold. Otherwise, there exists a partition n G II 
that contains at least three subsets of X. For every A G 7r, fix an arbitrary 
element xa G X \ A, define 11^ to be the set of the restrictions ^\au{x a } °f 
those partitions 7r' G II with the property that there exists some A' G 7r' 
with A U A' = X, and note that any such partition tt' can be recovered from 
its restriction 7rf , s x as it must consist of all subsets B in that restriction 
that do not contain xa and the complement of their union. Thus, it is not 
hard to see that, for every A G 7r, any two partitions of A U {xa} in Ha are 
compatible, that 1 + I^aI = 1-^1 holds, and that |v4U{xa}| < |^| holds 

for every A G n. Hence, by induction, 

|n| = i + ^2 l n ^l ^ 1 + X^ 3 ! 74 ! ~ 2 ) - 3n ~ 5> 

Ae-rr A£tt 

as required. ■ 

4 The algorithm for computing Cut*(D) 

In this section, we present our new algorithm for computing Cut*(D) called 
ComputeCutPoints(-D) which follows the recursive approach suggested by 
Lemma 13.11 This algorithm can be regarded as a speed-up of the algorithm 
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for computing outpoints presented in [7J, which, as mentioned in the intro- 
duction, also improves upon the run time of the algorithm presented in [2] . 
In Figure^ we present a pseudocode for this algorithm. Besides Cut*(D) the 
algorithm returns the set E^ and the auxiliary set A(Hd), which, for every 
split S = A\B G Hp, contains the 4-tuple (as, bs, D(as\B), D(bs\A)), where 
as G A and bs G B are fixed elements that are arbitrarily chosen during the 
course of the algorithm. 

To illustrate how our algorithm computes Cut*(D), consider the metric D 
presented in Figure [T^a). Suppose in Line 3 of the pseudocode in Figure El 
we select the element c. Consider the restriction D' of D to the subset 
X' := X \ {c}. A block realization of D' is presented in Figure IU(c). It is 
easy to check that the set of block splits of D' is 

s ' = {Wl{ & ; d > e }> WIK e}, {d}\{a, b, e}, {a, b}\{d, e}}. 

Note that the splits in E' are in one-to-one correspondence with the edges of 
the block realization in Figure [Tfc) . The set C := Cut*(D') consists of the 
Kuratowski maps in K(D') and one additional map / G R x ' with f(a) = 2, 
f(b) = 1, f(d) = 7 and /(e) = 3. Note that this map corresponds to the cut 
vertex v in Figure [Tfc) as f(x) equals the length of a shortest path from v 
to x in the block realization for every x G X' . 

Given C and E', the algorithm first computes the set E of block splits of 
D and the auxiliary set A (Lines 6-21). In our example it is easy to check 
that each of the splits in E' gives rise to precisely one split in E, that is, 

E = {{a}|{6, c, d, e}, {b}\{a, c, d, e}, {d}\{a, 6, c, e}, {a, b}\{c, d, e}}. 

Next the set C := Cut*(D) is computed (Lines 22-27) by first adding the 
maps Ja and fg for every S = A\B G E. For the metric D in Figure QJa), 
this yields, in addition to the Kuratowski maps k a , kb and kj, the 3 cutpoints 
(2,1,4,7,3), (3,2,3,6,2) and (8,7,2,1,3), where (x u x 2 , ■ ■ ■ , x s ) G M 5 rep- 
resents the map / G R x with (xi, x 2 , . . . , 25) = (/(a), fib), . . . , /(e)). Note 
that these cutpoints correspond to the 3 cut vertices in the block realization 
of D in Figure [T](b). For our example, the computation of C is completed by 
adding the Kuratowski maps k c and k e (Line 27). 

Theorem 4.1 Given a metric Dona set X with n elements, the algorithm 
ComputeCutPoints(-D) computes Cut*(D) in 0(n 3 ) time. 
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ComputeCutPoints(D) 



Input: a metric D on X 
Output: Cut*(D), E D , A(Z D ) 



1. if X = {x}, then return C := {k x }, E := and „4 := 0. 

2. Initialize C := 0, E := and .A := 0. 

3. Select iGl arbitrarily. 

4. Put X' := X \ {x}, and let D' denote the restriction of D to X'. 

5. Compute recursively C := Cut* (IT), E' := E D , and A' := A(E D /). 

6. for each 5" = A'\B' G S' do 

7. Put a s := og/ and 65 := 65/. 

8. Put A := A' U {x}, S := B' and extend S" to S := A\B. 

9. Compute D{a s \B) := D(a S /|B'). 

10. Compute := minlDfos/IA'), ^min (bsx + o 5 a — ax)}. 

A aeA 

11. if 5 is a block split of X, then 

12. Insert 5 into E and (a s , b s , D(a s \B), D(b s \A)) into A. 

13. Put A := A', B := B'U {x} and extend S' to S := A\B. 

14. Compute D(a s \B) := min-fD^/l -£?'), |min (asx + asb — ax)}. 

15. Compute D(b s \A) := D(b s >\A). 

16. if S is a block split of X, then 

17. Insert 5 into E and (as, bs, D(as\B), D(bs\A)) into A. 

18. Put 5 = {x}|X', a s := x and select b s G X' arbitrarily. 

19. Compute D(a s \X') and D(b s \{x}). 

20. if S is a block split of X, then 

21. Insert S 1 into E and (as, 65, Z?(as|X'), D(b s \{x})) into A 

22. for each S = A\B G E do 

23. Insert /a and f B into C. 

24. for each /' G C" do 

25. Extend /' to / G 1R X putting /(x) := max{xy - /'(y) : y G X'}. 

26. if / is a cutpoint of D, then insert / into C . 

27. for each x E X do insert fc x into C. 

28. return C, E and A. 

Figure 2: Pseudocode for our algorithm for computing Cut*(D). 
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Proof: We first show that our algorithm is correct. To do this we use in- 
duction on the size n of X. Our induction hypothesis is that our algorithm 
computes Cut*(D) and the set E^ of block splits of X correctly. If |X| = 1, 
there is nothing to prove. Now suppose that |Jf| > 1 holds. Let x be the 
element in X selected by our algorithm (Line 3), put X' :— X \ {x}, and 
let D' denote the restriction of D to X' (Line 4). By Lemma [3.1( i). the set 
Exj of block splits of X can be computed from the set E^>' of block splits 
of X'. By induction, the recursive call (Line 5) will correctly compute 
and, therefore, our algorithm will correctly compute E^ (Lines 6-21). Simi- 
larly, by Corollary 12.31 and Lemma I3~l7 ii) . the set Cut*(D) can be computed 
from Hjj and Cut*(D'). We have argued already that the computation of E# 
is correct and, again by induction, the recursive call (Line 5) will correctly 
compute Cut*(D'). Hence, our algorithm will correctly compute Cut*(D) 
(Lines 22-27). 

We next show that our algorithm has run time 0(n 3 ). We claim that 
an upper bound T(n) on the run time will satisfy the recurrence T(n) < 
T{n — 1) + 0{n 2 ). Using standard techniques for solving recurrences (see e.g. 
[1]), this yields T{n) G 0(n 3 ). So, it remains to show that all operations 
except those performed in the recursive call (Line 5) can be done in 0(n 2 ) 
time. 

We first focus on the computation of E^ from E^/ (Lines 6-21). Let 
S' = A'\B' be an arbitrary split in E^/. We can assume that D{as'\B') and 
D(b S '\A') are available from the 4-tuple (a s >, b s >, D(a s >\B'), D(b s >\A')) G A'. 
We want to check whether the split S = A\B = A' U {xjl-B' is a block 
split of X (Line 11). By Lemma 12.21 it suffices to check whether as > 
and asb + abs = asbs + ab holds for all a G A, b G B, using as = as> and 
bs = bs>. Note that, since 5" is a block split of X', it suffices to check whether 
asb + xbs = asbs + xb holds for all b G B, which can be done in 0(ri) time. 
Moreover, since D(as\B) = D(as>\B') and 

D(b s \A) = mm{D(b s >\A'),^mm{b s x + b s a - ax : a G A' U {x}}} 

hold (Lines 9-10), we can also compute as = D(as\B) + D{b s \A) — asbs in 
0(ri) time. 

To summarize, whether S is a block split of X or not can be checked 
in 0(n) time. Using completely similar arguments, it can also be shown 
that checking whether A'\B' U {x} is a block split of X (Line 16) can be 
done in 0(n) time. Note that, by Lemma 13.2^ there are 0(n) block splits 
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of D'. Thus, our algorithm will perform 0(n) iterations of the loop in Line 
6 and each iteration is completed in 0(n) time, yielding 0(n 2 ) in total for 
Lines 6-17. 

To finish the computation of we need to check whether the split 
S = {x}\X' is a block split of X (Lines 18-21). To do this, we fix as = x and 
choose an arbitrary 65 G X' . Then, we compute D(as\X') and D(bs\{x}), 
which can be done in 0(n 2 ) time, and check whether as = D(cts\X') + 
D(bs\{x}) — asbs > holds. We also check whether asb + xbs = asbs + xb 
holds for all b G X', which can be done in 0(n) time. This finishes the 
analysis of the time needed to compute £/> 

Next, we focus on the computation of Cut*(D) (Lines 22-27). We use 
a data structure Die to store the elements in Cut*{D) computed so far. 
Since, by Lemma [3721 \Cut*(D)\ G 0(n), the data structure Die can be 
implemented in such a way that inserting a single element of Cut*(D) into 
Die and, later on, checking whether an element of Cut*(D) has already been 
stored in Die both takes 0(n) time, see e.g. [llj . Moreover, we assume that, 
for every /' G Cut*(D'), the connected components of the graph Tfi have 
been computed and the cliques among them have been marked. 

So, first consider an arbitrary block split S = A\B G £d. If we have A = 
{x} and B = X', then we compute /y along with the connected components 
of Tf Y , marking the cliques among them, in 0(n 2 ) time for all Y G {A, B}. 
Next we consider the case that there exists some 5" = A'\B' G T, D > such that 
A = A' U {x} and B = B' (the following argument is completely analogous 
if A = A' and B = B' U {x}). Let a s G A' and b s G B' be the elements that 
we fixed for S in the course of the algorithm and let /a> and f B > be the maps 
in Cut*(D') associated with the split S'. Then we have 

f B (a) = D{a\B) = D{a s \B) - a s b s + ab s = D(a s \B') - a s b s + ab s = /b'(o) 

for all a G A' and 

f B (b) = D(b\A) -a s = a s b - D(a s \B) = a s b - D(a s \B') = f B ,(b) 

for all b G B — B', since D(as\B) = D(as\B') clearly holds. Hence, com- 
puting fs, the connected components of Tf B and marking the cliques among 
them can be done in 0(n) time, based on f B > and the connected compo- 
nents of Tf , . Similarly, if D(bs\A) = D(bs\A') holds, /a, the connected 
components of Y f A and the cliques among them can be computed in 0(n) 
time. Otherwise, that is, if D(bg\A) < D(bg\A') holds, the graph induced 
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by Tf A on X' is the disjoint union of two cliques with vertex sets A' and B', 
respectively. To see this, note that /a (a) + /a(o') > /a'(o) + fA'(a') > aa', 

+ = 2a 5 + + a 5 6' - 2£>(a 5 | J B) 

> a 5 6 + a s b' - 2D(a s \B>) = f B ,{b) + f B ,(b') > W, 

and f A (a) + f A (b) = ab = f A >(a) + f A >(b) holds for all a, a' G A' and b, b' G 
But then, also in this case, the connected components of Tf A and the cliques 
among them can easily be computed in 0(n) time. 

It remains to consider an arbitrary /' G Cut*(D') (Lines 24-26). Extend- 
ing /' to / (Line 25), that is, computing f(x) can be done in 0(n) time. 
Recall that we assume that the connected components of the graph Tfi and 
the cliques among them have been computed. From this information, we 
can compute in 0(n) time the connected components of Tf and determine 
which of them are cliques. Hence the loop in Line 24 will take 0(n 2 ) time, 
as required. Similarly, the loop in Line 27 will also take 0(n 2 ) time. This 
finishes the analysis of the run time of our algorithm and thus the proof of 
the theorem. ■ 
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